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ABSTRACT 

Recent X-ray observations of clusters of galaxies have shown that the entropy of the intracluster 
medium (ICM), even at radii as large as half the virial radius, is higher than that expected from grav- 
itational processes alone. This is thought to be the result of nongravitational processes influencing the 
physical state of the ICM. In this paper, we investigate whether heating by central AGN can explain the 
distribution of excess entropy as a function of radius. The AGN are assumed to inject buoyant bubbles 
into the ICM, which heat the ambient medium by doing pdV work as they rise and expand. Several 
authors have suggested that this "effervescent heating" mechanism could allow the central regions of 
clusters to avoid the "cooling catastrophe". Here we study the effect of effervescent heating at large 
radii. We find that the results are mainly sensitive to the total energy injected into the cluster. Our 
calculations show that such a heating mechanism is able to solve the entropy problem, provided that 
the total energy injected by AGN is roughly proportional to the cluster mass. The inferred correlation 
is consistent with a linear relation between the mass of the central black hole(s) and the mass of the 
cluster, which is reminiscent of the Magorrian relation between the black hole and bulge mass. 

Subject headings: galaxies: clusters: general — galaxies: active — dark matter — X-rays: galaxies: 
clusters — intergalactic medium — cooling flows 



1. INTRODUCTION 

Recent X-ray observations of both rich and poor clusters 
of galaxies have shown that there are problems in under- 
standing the temperature, gas density or, equivalently, the 
entropy profiles at large radii O.Sr^ir, where r^;, is the 
virial radius) (Ponman et al. 2003). The observed entropy 
at O.lraoo and at rsoo (where rjoo and r^oo correspond to 
radii within which the average overdensity is 200 and 500, 
respectively) is higher than that estimated from the purely 
gravitational interaction of gas with dark matter (Ponman 
et al. 2003 and references therein; see also Roychowdhury 
& Nath 2003). Earlier X-ray observations by Lloyd-Davies 
et al. (2000) had shown that entropy at O.lrsoo reached a 
"floor" for poor clusters and groups. However, recent re- 
sults of Ponman et al. (2003) have shown that the observed 
entropy is higher than the gravitational expectations for 
all clusters with emission- weighted temperatures (T) in the 
range 1-10 keV. The entropy at the much larger radius of 
rjoo is also found to be higher than expected from purely 
gravitational processes. 

Many theoretical models have been proposed to explain 
this phenomenon, including models that involve heat in- 
put from supernovae (Valageas & Silk 1999; Wu, Nulsen 
& Fabian 2000), gas cooling (Bryan 2000; Voit & Bryan 
2001; Muanwong et al. 2002; Wu & Xue 2002; Dave, Katz 
& Weinberg 2002, Tornatore et al. 2003), accretion shocks 
(Tozzi & Norman 2001; Babul et al. 2002) and quasar out- 
flows (Nath & Roychowdhury 2002). 

Observations have also revealed the presence of X-ray 
deficient bubbles in the inner regions of many cooling flow 
clusters, e.g., the Hydra A cluster (McNamara et al. 2000), 
Abell 2052 (Blanton et al. 2001, 2003), Abell 2597 (McNa- 



mara et al.2001), AbeU 4059 (Heinz et al.2002), Abell 
2199 (Johnstone et al.2002), and others. These bub- 
bles are characterized by low X-ray emissivity, implying 
low density compared to the ambient medium. In most 
of these cases, the cavities are clearly coincident with 
the radio lobes of the AGN in the cluster center. How- 
ever, some clusters also exhibit cavities with weak or un- 
detectable radio emission (known as "ghost bubbles" or 
"ghost cavities") located far away from the cluster cen- 
ters, like Perseus (Fabian et al.2000), MKW3s (Maz- 
zotta et al. 2002), and Abell 2597 (McNamara et al. 2001). 
These bubbles are also believed to be filled with relativistic 
plasma or buoyant gas deposited by jets from the central 
AGN, and are thought to rise through the ICM subson- 
ically due to buoyancy. As they rise, radiative and adi- 
abatic losses reduce the energy of the relativistic plasma 
inside the bubbles, resulting in a very low radio flux. The 
discovery of these bubbles and their detailed observational 
study have stimulated theoretical studies of the impact of 
these bubbles on the intracluster medium. 

The evolution of these bubbles has been studied exten- 
sively in connection with the cooling flow catastrophe in 
the centers of clusters. It has been found recently that 
simple cooling flow model of clusters, which predict that 
the temperature of gas in the central regions of clusters 
should be very low (much less than 1 keV), are mostly in 
conflict with the observed temperature profiles (Peterson 
et al. 2001; Allen et al. 2001 and many others). Several au- 
thors have addressed whether the absence of very cold gas 
(below ~ 1 keV) can be explained by AGN energy input via 
outflows and bubbles (see, e.g.. Tabor & Binney 1993, Bin- 
ney & Tabor 1995, Ruszkowski & Begelman 2002, Chura- 
zov et al. 2002, Kaiser & Binney 2003, and Briiggen 2003a 
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for semi-analytic models; and Quilis ct al. 2001, Reynolds 
et al. 2002, Briiggcn 2003b, and Ruszkowski, Briiggen & 
Begelman 2004, Basson & Alexander 2003, Omma et al. 
2003a, Omma & Binney 2003b for numerical simulations). 

In this paper, wc explore the possibility of heating the 
intracluster gas at large radii via the "effervescent heat- 
ing" mechanism (Begelman 2001, Ruszkowski & Begelman 
2002). Since the excess entropy requirements at different 
radii are different, the questions we set out to answer are 
the following: 

1. Is it possible to satisfy simultaneously the entropy 
observations at two fiducial radii (O.lrjoo and r^oo) 
with a single central heating source for all clusters? 
We explore the parameter space of time-averaged 
jet luminosity (L) and the time for which this ac- 
tivity continues, theat- 

2. How does the entropy profile of the cluster gas 
evolve with time when it is being heated and cooled 
simultaneously? 

The paper is organized as follows. In § 2 we describe 
our model, including initial conditions, details of the effer- 
vescent heating mechanism, our prescriptions for cooling 
and convection, and the simulation method. The results 
and discussion are presented in §§ 3 and 4, respectively. 
We discuss caveats of our model in § 5 and summarize our 
conclusions in § 6. 

2. MODEL 

Our model of the ICM assumes that the thermal gas re- 
mains in hydrostatic equilibrium as it is heated by buoyant 
radio bubbles, originating from the central AGN, via the 
"effervescent heating" mechanism. Here, we begin with a 
description of the background dark matter potential of the 
cluster. 

2.1. Dark matter density profile 

We assume that the dark matter density profile of the 
cluster, pa,„{r), is described by a self-similar form as sug- 
gested by many high resolution AT-body simulations. The 
profile is expressed in terms of a characteristic radius, r^, 
by 

(see, e.g., Komatsu & Seljak 2002), where is a charac- 
teristic density and x = r/r^. We define a dimensionless 

"concentration parameter" c = r^^Jr^, where the virial ra- 
dius r^i, is calculated from the spherical collapse model 
(Peebles 1980) assimiing overdensity Ac{z = 0) = 100 
(Komatsu & Seljak 2002). The characteristic density pe is 
then given by 



47rr^. 



where M^;, is the cluster virial mass, 
parameter c can be approximated by 



The concentration 



c= gf 



-0.13 



according to numerical simulations by Bullock et al. (2001) 
(we use h = 0.71). 



The above set of equations specifies the dark matter 
density profile of a particular cluster as a function of its 
virial mass. Next, we turn our attention to the density 
profile of the gas in hydrostatic equilibrium with this dark 
matter distribution. 

To compare our results with observations, we present 
our results in terms of the radii rsoo and r^oo, where 
the overdensities are 200 and 500, respectively. Both 
radii are functions of the cluster mass. The ratio 
r-Ma/r-iaa ~ 0.65 — 0.67 for the range of masses we have 
considered. 



2.2. Initial configuration of intracluster gas 

We use the "universal temperature profile" (Loken et 
al. 2002) as the initial temperature profile of gas in hy- 
drostatic equilibrium to determine the initial density pro- 
file. This profile is different from the commonly used 
self-similar profile, which assumes that the intracluster 
gas density distribution scales as the dark matter den- 
sity profile with a constant of proportionality A = ^ 
(i.e., Pgaa = ./bAim)- We use the "universal temperature 
profile" instead of the self-similar profile because the for- 
mer is claimed to be the result of gravitational processes 
alone. Moreover, it does not have the unrealistic turnover 
of the temperature profile in the inner regions of the clus- 
ter, which one finds in the self-similar profile (Roychowd- 
hury & Nath 2003). 

The initial temperature profile (normalized by the 
emission- weighted temperature (T)), in the radial range 
0.04rvir <r< rvir, is given by 



r„(r) _ b 
"(T)~ ~ {l + r/af 



(4) 



where h = 1.33, a = r^i,/1.5, and S = 1.6. To determine 

the emission- weighted temperature from the cluster mass, 
we use a relation that arises from adiabatic evolution of 
the gas in the cluster. Afshordi & Cen (2002) have shown 
that the observations by Finoguenov et al. (2001) of the 
M5oo-(T') relation in clusters can be understood from grav- 
itational processes alone. We therefore use this empirical 
relation {M^oo-{T)) derived by Finoguenov et al. (2001): 



M.nn = 2.64 X 10 



13 



/fcb(T)y 

VlkeV/ 



78 



(5) 



where fct is the Boltzmann constant and M^^o has been 
calculated self-consistcntly by taking the total mass within 
the radius where the over-density is (5 > 500. 

The equation of hydrostatic equilibrium for gas in a clus- 
ter with temperature T<,(r) and density Pgaso(?') is 



Pga.o(r) dr 



,-Pgaso('')) 



GM(< r) 



(6) 



where Pg^so = {Ps^so/lJ'm^)KTo, M{< r) is the total mass 
inside radius r, and p, and rUp denote the mean molecular 
weight (/i = 0.59) and the proton mass, respectively. The 
boundary condition imposed on the solution is that the 
gas-fraction, /g^, = Mgas/M^m = 0.13333, within raoo is 
universal, as recently found by Ettori (2003) for a sample 
of low- and high-redshift clusters. 
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2.3. Heating, cooling and convection 
2.3.1. Heating 

In the "effervescent heating" model, the central AGN is 

assumed to inflate buoyant bubbles of rclativistic plasma 
in the ICM (Begelman 2001, Ruszkowski & Begelman 
2002). The timescale for the bubbles to cross the clus- 
ter (of order the free-fall time) is shorter than the cool- 
ing timescale. Since the number flux of bubbles is largo, 
the flux of bubble energy through the ICM approaches a 
steady state. This implies that details of the energy in- 
jection process such as the number flux of bubbles, the 
bubble radius or size, filling factor and rate of rise do not 
affect the average heating rate. 

It is assumed that the rclativistic gas docs not mix with 
the ICM very efficiently. Under such conditions the bub- 
bles can expand and do pdV work on the ambient medium 
as they rise in the cluster pressure gradient. We assume 
that this work converts internal energy of the bubbles to 
thermal energy of the intracluster gas within a pressure 
scale height of where it is generated. Then, in a steady 
state (and assuming spherical symmetry), the energy flux 
carried by the bubbles varies as 



n(r) 



(7b-l)/7b 



(7) 



where Pb{r) is the partial pressure of buoyant gas inside 
the bubbles at radius r and 7^ is the adiabatic index of 

buoyant gas, which wc have taken to be 4/3. This formula 
includes the effects of adiabatic bubble inflation. Assum- 
ing that the partial pressure inside these bubbles scales 
as the thermal pressure of the ICM, the volume heating 
function H can be expressed as 



n ~ -r2/i(r)V • (fF,) 



= -h{r)( 



V R) / r d\nr 



(8) 



where Pg is some reference pressure and h{r) is the nor- 
malization function 



h{r) 



(1 - exp(-r/ro)) q 



(9) 



In equation (9), {L)is the time- averaged energy injection 
rate and rg is the inner heating cut-off radius which is 
determined by the size of the central heating source. In 
our calculations To is taken to be O.Olrvi,, which is ~ 15 - 
20 kpc depending on the cluster mass M^i{= Myir). The 
normalization factor q is defined by 

I min " 



where rn 



^ r2oo- 



2.3.2. Cooling and convection 

To calculate the volume cooling rate, wc use a fit to 
the normalized cooling function A^ (T) for a motallicity of 
Z/Zq = 0.3, as calculated by Sutherland & Dopita (1993). 
This cooling function incorporates the effects of free-free 
emission and line cooling. The fit is borrowed from Nath 



(2003). Thus, the volume cooling rate is T = n'^A^{T), 
where = 0.875(p/mp) is the electron density. 

The convective fiux i^conv is given by the mixing length 
theory. 



otherwise, 

(11) 

where g is the gravitational acceleration, /„ is the mix- 
ing length, s = (7 - ly^K/inrUp) ln{P^^Jp2J is the gas 
entropy per unit mass, and Cp = "fk-^/Kj — l)/iTOp] is the 
specific heat per unit mass at constant pressure. We use 
— min[0.3Pg,i,/(/9gasS'), r], where r is the distance from 
the cluster center. 

2.3.3. Evolution of the intracluster gas 

As noted earlier, the gas is assumed to be in quasi- 
hydrostatic equilibrium at all times since the cooling is 
not precipitous at these radii and the heating is mild. The 
gas entropy per particle is 



S = const + 



1 



-fcbln(cr). 



where a = P^^^l Pg^s the "entropy index" and 7 is the 
adiabatic index. The particle number density of the gas, 
n, is given by n = p^J/ nm^. 

During each timestep Ai, the entropy of a given mass 
shell changes by an amount 



AS' = 



1 



Act 



1 



(W - r - V • F,„„OAt. (13) 



7 — 1 "a uT 
Incorporating the expressions for heating, cooling, and 
convection, the entropy increment for each mass shell for 
a timestep At becomes 



Aa{M) 



3 P.. 



-At 



1 (ilnP„ 



r d In r 



1 . 2„ 

(r F 

ar 



(14) 



where n„ and Pg,,^ are the current electron number density 
and pressure of the ICM, respectively. 

Thus, the entropy index of each mass shell of gas due to 
heating and cooling after a time At becomes 



a_(M) = ao(M) + Aa(M) 



(15) 



where (To(Af) = Pga.so//5g„o the default entropy index. 
The system relaxes to a new state of hydrostatic equilib- 
rium with a new density and temperature profile. After 
updating the function a{M) for each mass shell, we solve 
the equations 



rfPga, _ GM{< r) 
dM 



47rr* 



dr _ 1 / Pga3(M) N 
dM 47rr=\ a{M) J 



(1/7) 



(16) 



(17) 



to determine the new density and temperature profiles at 
time t + At. The boundary conditions imposed on these 
equations are that (1) the pressure at the boundary of the 
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cluster, Tout, is constant and is equal to the initial pres- 
sure at Taoo, i-e., P(ro„t) = constant = Pga=o(''2oo), and 
(2) the gas mass within Tout at all times is the mass con- 
tained within Tsoo for the default profile at the initial time, 
i.e., A4(r„t) = A4o(?'2oo) = 0.1333Afdm(r2oo)- It is impor- 
tant to note here that rout increases as the cluster gas gets 
heated and spreads out. 

The observed gas entropy <S(r) at O.lraoo and at r^oo is 
then calculated using 

Sir)^Tir)/n^J'ir). (18) 

The updated values of a{M) and pressure of the ICM 
PgB.s{r) are used to calculate the heating and cooling rates 
and the convcctivc flux for the next time step. This is 
continued for a duration of ^hoat- After that the heating 
source is switched off, putting H = 0. The cooling rate 
and convective flux continue to be calculated to update 
the function a{M) at subsequent timesteps, and the hy- 
drostatic structure is correspondingly evolved for a dura- 
tion of tu — theat, where ta is the Hubble time. Note that 
Tout decreases during this time since the intracluster gas 
loses entropy and shrinks. The only free parameters in 
our calculation are the energy injection rate (L) and the 
time fi„,„, over which the "effervescent heating" of the ICM 
takes place. After evolving the gas for the total available 
time, ^ 1.35 x 10"'^° years, we check whether the entropy 
at O.lrsoo and r^oo matches the observed values, and ad- 
just parameters accordingly. In this way, we explore the 
parameter space of (L) and for different cluster masses 
so that the entropy (after 1.35x10^" years) at O.lraoo and 
Tsoo matches the observed values. 

For numerical stability of the code, the convection term 
is integrated using timesteps that satisfy the appropriate 
Courant condition. The Courant condition for convection 
is 

< l^'^Wj^- (19) 

The timesteps, At, used in equation (14) to update the 
entropy of the gas and calculate its pressure, temperature 
and density profiles always obey the above Courant con- 
dition (Ruszkowski & Begelman 2002; Stone, Pringle & 
Begelman 1999). 

The calculations presented in this paper were also per- 
formed using the fully time-dependent ZEUS code (Clarke, 
Norman & Fiedler 1994) and the results obtained were con- 
sistent with the ones presented here. However, we decided 
to use our quasi-hydrostatic Lagrangian code as it allowed 
us to search the parameter space more efficiently. 

3. RESULTS 

In this section, we discuss our results for cluster evo- 
lution due to heating, cooling and convection. The gas 

is heated for a time ihcat and cooled simultaneously. Af- 
ter this time, the heating source is switched off. The gas 
is then allowed to cool radiatively until a total simulation 

time of tH=l-35x IQ-'^'' years has elapsed. The final entropy 
values at O.lraoo and r^oo are compared with the observed 
ones. 

In Figure 1, the observed entropy values and their er- 
rors at O.lraoo and at r^oo are plotted as a function of the 



emission-weighted temperature, (T) , of the cluster (Pon- 
man et al. 2003). We have done a best-fit analysis on these 
data points to estimate the entropy requirements for the 
sample of clusters of masses ranging from 10^^ to 2 x 10^^ 
Mq. The results of our analysis arc the shaded regions 
in the two panels of Figure 1. The solid line through the 
center of the shaded region is the best-fit curve with the 
lines bounding the shaded region being the 1 — a errors on 
the best-fit entropy values. 

Figure 2 shows the time evolution of scaled entropy pro- 
files of a cluster of mass M^i = 6 x 10^'' Mg with and 
without convection for (L) = 3x10^^ erg s~^. We use the 
same method of emissivity weighting as in Roychowdhury 
& Nath (2003) to calculate the average quantities. The 
entropy profiles are plotted in time-steps of 5 x 10^ years. 
They are seen to rise with time as the ICM is heated. 
Then, after the heating is switched off (after theat = 5 x 10^ 
years), the gas loses entropy due to cooling and the pro- 
files are seen to fall progressively until the Hubble time 
is reached. The inclusion of convection flattens the nega- 
tive gradient in the scaled entropy profiles in the central 
regions of the cluster (within 0.2r2oo, as seen in the dif- 
ference of the entropy profiles between the left and right 
panels in Figure 2. The plot that includes convection (left 
panel in Figure 2) shows that the gas develops a flat en- 
tropy core in its central regions after the heating source 
has been switched off and the gas has cooled. 

Figure 3 shows the evolution of density and tempera- 
ture of the ICM for a cluster of mass 6 x 10^^ Mq and 
for a luminosity of (L) = 3 x 10^^ erg s~-^. The gas den- 
sity decreases with time during the heating epoch, and 
increases due to radiative cooling and convection after the 
heating source is switched off. It is interesting to note 
that the changes in density are minimal beyond 0.2r2ooi 
and that convection plays an important role in regulating 
the density profiles after the heating source is switched off 
for radii r < 0.2r2oo- We note that for clusters with lower 
emission- weighted temperatures, the effects of heating and 
convection are seen at larger radii. 

We now discuss the permitted range in the total energy 
injected into the cluster, = (L) x theat, required to 

match the observed entropy as a function of the cluster 
mass. Figure 4 shows the spread in E^g„ as a function of 
the cluster mass for two different values of theat- The re- 
gion bounded by thin solid lines is the permitted range in 
energy for fi„,at = 5 x 10^ years. This region includes an 
area shaded with solid vertical lines, which corresponds to 
the values of E^g„ that satisfy the entropy requirement at 
O.lrjooj and another area shaded with solid oblique lines 
which satisfies the entropy requirement at Tsoo. Similarly, 
the region bounded by thin dotted lines in Figure 4a shows 
the permitted spread in energy for theat = 5 x 10^ years. 
This region includes an area shaded with dots which cor- 
responds to the values of E.^^^ that satisfy the entropy 
requirement at 0.1r2oo and another area shaded with long- 
dashed horizontal lines, which satisfies the entropy require- 
ment at Tsoo- For cluster masses above 9 x lO^^M© there is 
no lower limit on the injected energy from the entropy mea- 
surements at rsoo- The region corresponding to this situ- 
ation is marked by by horizontal dashed lines and oblique 
lines for cluster masses above 9 x lO^^M©. 

Figure 5 shows the permitted total injected energy range 
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as a function of the mass of cluster for heating times be- 
tween theat = 5 X 10^ years and theat = 5 x 10^ years. 
Here the entropy is required to match observations at both 
O.lrsoo and r^oo- The thick sohd Une represents a hnear 
relation between the total energy injected to the cluster 
by AGN and the mass of the cluster (see next section for 
more details). 

4. DISCUSSION 

As is clear from Figures 4 and 5, it is possible to heat 
the ICM with a single central AGN to match the entropy 
requirements at both O.lrjoo and r^oo- However, in order 
to match the entropy at both radii, the total injected en- 
ergy -Eagn, for a given value of th„,jt <C t^, must be tightly 
constrained. In fact, our calculations have shown that for 
any value of theat < in, i-e., for any heating time (or AGN 
lifetime), it is always possible to satisfy the entropy ob- 
servations at both radii with a single value of the luminos- 
ity (L) . This is different from the cooling flow problem, 
where (L) must be finely tuned to match the cooling rate 
(Ruszkowski & Begelman 2002), because cooling effects on 
large scales are rather mild and, thus, the results depend 
mostly on the total injected energy E^^^ = (L) x fheat- 
Thus, if we can fit the observed entropy values for just 
one pair (L) and theat, we can do so for a wide range of 
such pairs. Another manifestation of this fact is that the 
plots in Figure 4 arc similar in shape but arc just offset 
by a factor of a few. As the cooling effects are relatively 
mild, the modest differences arise because shorter heating 
times lead to higher temperatures and convective trans- 
port cannot "catch up" with the energy supply. As the 
entropy excess to be explained is known from observations 
and is fixed, more heat has to be injected for shorter heat- 
ing times. Nevertheless, the results are mostly sensitive to 
the total energy input from the black hole, rather than to 
(L) and ihcat separately. As a consequence, satisfactory fits 
can be obtained as long as the total injected energy falls 
within a relatively narrow range of values, which depends 
on the cluster mass (Figure 5). 

As an illustrative example of such a relation, the thick 
solid line plotted in Figure 5 shows a linear relation be- 
tween E^g^ and cluster mass, M^i, that is consistent with 
the requirements from entropy observations for a heating 
time of 5 X 10^ < theat < 5 x 10^ years. We point out 
that the effective heating time may be longer than the 
integrated AGN lifetimes. Heating of the cluster may oc- 
cur more gradually as heat gets distributed on a timescale 
very roughly proportional to the sound crossing time from 
the cluster center to a given radius. Thus, thoat should be 
interpreted as an upper limit on the AGN lifetime. 

Let us assume that the bulge mass of the central AGN 
is related by a factor / to the cluster mass as M,, ^ f M^i 
and that the mass of the black hole is related to the bulge 
mass of the AGN host by Mt,h ~ 1.5 x IQ-^M^, (Haring & 
Rix 2004; Fcrrarcsc & Merritt 2000, Gebhardt et al. 2000). 
This means that the mass of the black hole is related 
to the cluster mass by Mbh ~ 1-5 x 10~^/_2Afci, where 
/_2 = //O.Ol. It is generally thought that AGN are fu- 
elled by accretion onto the supermassive black hole. If 
we require that a fraction 77 of the mass of the black hole 
is converted into energy, then the relation between the 
black hole mass and the total injected energy becomes 



Mbh 2.8 x 10^f]-lEe2MQ, where t]o.2 = Jy/O-S. The 
relation between the total injected energy and the clus- 
ter mass is Eq2 5.4?7o.2/-2-Mi4. The linear relation in 
Figure 5 corresponds to /-2??o.2 = 0.25. 

Note that for longer heating times the luminosity con- 
straints become less stringent and lower black hole masses 
can explain the observed trends. Finally, as we consider 
heating on large scales, all galaxies in the cluster that go 
through an active phase will contribute to the overall en- 
ergy budget of the cluster gas. For example, substructure 
in the cluster could contain galaxies with sizable bulge 
components, each of which may contain a supermassive 
black hole. Therefore, the constraint on the black hole 
mass obtained above should be interpreted as a sum of all 
black hole masses that contribute to heating rather than 
the mass of an individual black hole. This could also lower 
the required efficiency of individual black holes that con- 
tribute to cluster heating. 

Finally, we note that cooling and convection play impor- 
tant roles in controlling the heating mechanism so that the 
entropy profiles broadly match the observed entropy pro- 
files in clusters (Ponman et al. 2003). Notably, in the later 
stages of evolution of the gas, after the heating source is 
switched off, convection flattens the entropy proflles in the 
central regions of the cluster. Moreover, an entropy core 
seems to develop in the final stages of the evolution of the 
ICM. This is consistent with the observed entropy profiles, 
which show cores at r < O.lrsoo (Ponman et al. 2003). 
Our entropy profiles do not show steep positive gradients 
as observed in the entropy profiles of two groups of galax- 
ies (Mushotzky et al. 2003). However, these groups have 
masses smaller than the ones considered here. Moreover, 
since the temperature profiles show a temperature gradient 
even though the entropy is flattened, thermal conduction 
might be able to conduct heat out from the central re- 
gions of the cluster and reduce temperatures there. This 
will happen at roughly constant pressure and, therefore, 
central densities will slightly increase. Thus conduction (if 
not fully suppressed) may help to remove the core in the 
center. Analyzing the effects of conduction is not within 
the scope of this paper. 

5. CAVEATS 

In this section we make explicit the caveats of our model. 
We note that our model assumes that all of the injected 
energy goes into effervescent heating. The amount of en- 
ergy that is available to effervescent heating, as opposed 
to the energy that is "lost" to bubble creation, can be 
estimated from the following simple analytic argument. 
The available effervescent energy equals the energy lost 
to bubble expansion as it rises buoyantly, fp PdV, where 
Po is the pressure at the transition region from the bub- 
ble formation region to the buoyant (effervescent) phase. 
To get the maximum energy available we set the upper 
limit of the integral to zero pressure. Assuming adiabatic 
evolution of the gas inside the bubbles (low density, low 
radiative losses) and mass conservation in the bubble one 
has dV = il/j){Po/P)^/'^VodP/P, where Vo is the volume 
of the bubbles at the transition region mentioned above. 
On integrating, effervescent energy is SPqVo for 7 = 4/3 
(appropriate to the bubble interior). Thus, since the en- 
ergy fraction lost to bubble creation at constant pressure 
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is just PoVo, approximately only 25% of the injected en- 
ergy goes to bubble creation, while 75% is available for 
effervescent heating. We note that if bubble creation oc- 
curs at constant pressure, then the injected energy is con- 
verted to pdV work on the ICM. The exact details of the 
dissipation of this energy are as yet unknown, but presum- 
ably this process occurs with efficiency lower than 100%. 
Thus, in the regime where bubble inflation occurs at con- 
stant pressure, the above estimate of energy lost to the 
bubble creation is probably a lower limit. However, if the 
bubble inflation in the cluster center occurs supersonically, 
the energy required to inflate the bubble goes into shock 
heating. Any shocked, high entropy gas is buoyant itself, 
and acts in a fashion similar to the case described by the 
effervescent heating function the difference being that a 
small fraction of the total mass from the innermost radial 
shells is removed. 

We note that our treatment is simplified as we neglect 
cosmological evolution. In the hierarchical picture, tem- 
perature increases with merging as small objects merge 
to form large objects with higher masses. Also, merging 
can only increase the entropy. Therefore, heating of the 
smaller units before they merge, and therefore when they 
have lower virial temperatures, should reduce the total 
energy requirement according to the second law of ther- 
modynamics {dE = TdS). As clusters observed today 
{z ~ 0) were most likely formed at z < 0.5 (Kitayama 
& Suto 1996, Balogh, Babid & Patton 1999) and AGN 
could switch on before the formation epoch of these clus- 
ters, lower energy injection could be required to explain 
observations at z ~ 0. 

We also note that the hydrodynamics equations have 
been integrated up to r2oo- It is likely that the buoyant 
energy transport would be inhibited by the cluster accre- 
tion shock at such large radii. However, excess entropy 
used to constrain our model is measured at smaller radii 
and we do not expect this to affect our results significantly. 

Our simulations assume that heating is instantaneous. 
This is an approximation as there should be a delay be- 
tween the onset of AGN activity and heating at a given 
radius. We have chosen the shortest heating duration to 
be at least the sound crossing time (or the dynamical time 
of rising bubbles) at the outer radius where we compare 
the computed and observed entropies (rsoo)- Thus, for the 
shortest heating time, our assumption at 0.1r2oo should be 
quite reasonable but at rsoo it would be less accurate but 
still physically conceivable. For longer heating times our 
results at both radii should not be affected by this approx- 
imation. 

Related is the issue of the constraints on the source lu- 
minosity. For shorter heating times, higher luminosities 
are required to assure the same energy injection. Thus, 
the required luminosities are more feasible for longer heat- 
ing times. We note also that if the heating is supplied 
by more than one AGN, then the luminosity requirements 
on an individual source would be lower. We also reiterate 
the point made earlier that if the effects of cosmological 
evolution had been fully taken into account and AGN had 
supplied heating at earlier epochs, then the energy require- 
ments presented here would have been reduced. Thus, the 



source luminosities would also be lower for a given heat- 
ing time. This is due to the fact that the required excess 
entropy at higher redshift would be reduced. 

6. CONCLUSIONS 

We have studied whether "effervescent heating" of clus- 
ter gas by a central AGN can resolve the entropy problem 
in clusters of galaxies. In this model, the AGN (or a group 
of AGNs in the central region) injects bubbles of buoyant 
gas, which heat the ICM. The mean volume heating rate 
due the bubbles is a function of the ambient pressure and 
a time-avcragcd energy injection rate to the ICM, but not 
of the detailed properties of the bubbles or their evolution. 
We have also included the effects of radiative cooling and 
convection. We assumed that heating continues for thea,t, 
the duration of heating, and have studied the resulting 
evolution of the gas assuming it to be in quasi-hydrostatic 
equilibrium for the Hubble time. The only free parameters 
of this model are {L) , the energy injection rate to the ICM, 
and fheaf The main results of our study are summarized 
as follows: 

1. We find that there are allowed values of (L) for 

which it is possible to match the entropy observa- 
tions of clusters even at large radii (both O.lrjoo and 
rgoo) with a single central AGN for a large range 
of iheat (5 X 10^ < theat < ycars) and cluster 
masses (5 x 10^^ - 2 x lO^^Mg). 

2. Convection plays an important role in removing 

negative entropy gradients produced by heating in 
the central regions of the cluster. 

3. The results are mostly sensitive to the total energy 
Eagn injected into the cluster by AGN (as cooling 
effects arc relatively mild). The model predicts that 
the total injected energy ii^agn ={L) xAt, required 
to satisfy observational entropy constraints, should 
be correlated with the cluster mass. This require- 
ment is consistent with a linear relation between 
the mass of the central black hole(s) and the mass 
of the cluster, which is reminiscent of the Magorrian 
relation between the mass of the black hole and the 
bulge mass of the host galaxy. 
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Temperature (keV) Temperature (keV) 
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Fig. 1. — Gas entropy (defined here as T/nJ ) as a function of emission-weighted temperature (T) , at radii 0.1r2oo (left panel) and rsoo 
(right panel). The data points are from Ponman ct al. (2003). The solid line in the center of the shaded region is the best-fit to the data 
points. The two solid lines bounding the shaded region are the 1 — ct errors on the best-fit values of entropy. The dashed line is the predicted 
entropy due to gravitational interactions alone (from Roychowdhury & Nath 2003). 
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Fig. 2. — Scaled entropy profiles as a function of scaled radius for a cluster of mass 6 X 10^'' Mq heated by an AGN with (L) = 3 X 10'*''' erg 
s""*^, with convection (left panel) and without convection (right panel). The scaled entropy profiles are plotted at intervals of 5x10^ years. 
They are seen to rise as the gas is heated and then fall as the gas cools. In both cases theat = 5 x 10^ years. Initial states correspond to the 
lowest curves. 
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Fig. 3. — Gas density (left panel) and temperature (right panel) profiles as a function of scaled radius (r/r2oo), for the cluster model shown 
in the left panel of Figure 2. It is seen here that radiative cooling lowers the temperature, thus increasing the density after the heating source 
is switched off. Initial density and temperature profiles correspond to the highest and lowest curves, respectively. 
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Fk;. 1. — Allowed ranges of total injected energy -Eagn = {L) X thoat required to match the observed entropy at 0.1r2oo and/or rsoo, for 
two different values of (heat , as a function of cluster mass. The width of each region corresponds to the 1 — cr errors plotted in Figure 1. The 
region bounded by thin solid lines corresponds to the allowed range in energies for thcat = 5 X 10** years. Within this bounded region, the 
area shaded with vertical solid lines is the spread in energy that satisfies the entropy requirement at 0.1r2oo and the area shaded with oblique 
solid lines is the spread in i?agn that satisfies the entropy requirement at rsoo. Only the overlap region satisfies the observations at both radii. 
The region bounded by thin dotted lines corresponds to the allowed range in energies for thoat = 5 X 10^ years, in this region, the dotted area 
is the spread in energy that satisfies the entropy requirement at 0.1r2oo and the area shaded with horizontal long-dashed lines is the spread 
in £agn that satisfies the entropy requirement at r^oo- 
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100 




Fig. 5. — This figure shows the permitted total injected energy range as a function of the cluster mass for ICM heating times between 
*heat = 5 X lO^yr (upper envelope) and theat = 5 x lO^yr (lower envelope). The shaded region corresponds to values of Sagn that are able 
to match the entropy observations at both 0.1r2oo and rsoo- The thick solid line represents a linear relation between the total energy injected 
into the cluster by AGN and the mass of the cluster. The permitted parameter space comes from the sum of permitted regions in Figure that 
satisfy the entropy constraints at both radii for fixed theat ■ 



